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Abstract 

The contributions from the color-octet components \(cb)s( 1 So)g) and \(cb)g( s Si)g) to the /is c 
or Xb c (the P-wave B c meson) hadroproduction are estimated in terms of the complete 0(a 4 ) 
calculation. As necessary inputs in the estimate, we take the values of the octet matrix elements 
according to the NRQCD scaling rules, and as a result, we have found that the contributions to the 
P-wave production may be the same in order of magnitude as those from the color-singlet ones, 
|(c6)i( 1 Pi)) and |(c6)i( 3 Pj)) (J = 1,2,3). Especially, our result indicates that the observation of 
the color-octet contributions to the P-wave production in the low transverse momentum region is 
not very far from the present experimental capability at Tevatron and LHC. 
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I. INTRODUCTION 

B c meson has been observed experimentally , and within theoretical uncertainties 
and experimental errors the observations are consistent with the theoretical predictions. 
In view of the prospects in B c physics at the Fermilab Tevatron II and LHC, B c physics 
is attracting more and more attentions. The hadronic production of the B r meson (its 
ground state mainly), has been studied quite well Especially, a 

generator BCVEGPY for hadronic production of B c meson has been available, which 

can be easily complimented to the PYTHIA environment and enhances the efficiency 
to generate the full B c events greatly in comparison with using PYTHIA itself. 

In the framework of effective theory of NRQCD jl^l . I3I ] , a heavy quarkonium is consid- 
ered as an expansion of various Fock states. The relative importance among those infinite 
ingredients is evaluated by the velocity scaling rule. A similar idea can be applied to the 
double heavy meson system too. That is, the physical state of B c , B*, h Bc and Xb c mesons 
can be decomposed into a series of Fock states: 

\B C ) = O^KcbUCSo)) + 0(v 1 )\(cb) 8 CP 1 )g) + 0(v 1 )\(d>) 8 ( 3 S 1 )g) + ■■■ 

\B*) = O(v )\(cb) 1 ( 3 S 1 )) + 0(v l )\(cb) 8 CPj)g) + 0(v 1 )\(d>) 8 ( 3 S l )g) + ■■■ (1) 

and 

\h Bc ) = O(i; )|(c6) 1 ( 1 P 1 )) + 0(v 1 )\(cbU 1 So)g) + 0{v l )\{ci)^Pj)g) + ■■■ 
Ixl) = Oiv^KcbU'Pj)} + 0(v 1 )\(cb) 8 ( 3 S l )g) + O^KcbU'Pjg) + ■ • • , (2) 

here v is the relative velocity. Note that throughout this paper we use the symbols hs c 
and Xb c to denote the four physical P-wave states {hs c denotes the P-wave state with 
the dominant color-singlet state (cfo)i( 1 Pi), while Xb c denotes the P-wave states with the 
dominant color-singlet states (c6)i( 3 Pj) as described in Eq.(J2J)) instead of B* JL=l that is 
used in Ref. (14|. Here the thickened subscripts of the (cb) denote for color indices, 1 for 
color singlet and 8 for color-octet; the relevant angular momentum quantum numbers are 
shown in the parentheses accordingly. According to the velocity scaling rule, the probability 
of each Fock state in the expansion is proportional to the scales in a definite power of v, 
squared the indicated one as the above. Generally, the leading Fock states of the B c and 
B* states are |(c6)i( 1 S'o)) and |(cfe)i( 3 S'i)) respectively, whose probability in respect of the 



relative velocity are set to be the order of 0(v°). Also for the production of the B c and B* 
states, there is no problem at all, one can be sure that the contributions from leading Fock 
state (s) are dominant. 

As implied in equations Eq.(j2J), at leading and next leading order, the production of P- 
wave B c mesons may involve in the Fock states: |(c6)i( 1 Pi)), \(cb) 8 ( 1 S )g) , |(c6) 8 ( 3 F/)g) ■ ■ •; 
|(c6)i( 3 Fj)), |(c6) 8 ( 3 Si)(7), |(c6) 8 ( 1 P 1 )(7) • • •. For color-singlet components (P-wave), the long 
distance non-perturbative matrix elements can be directly connected with the first derivative 
of the wave functions at the origin, which can be computed via potential models ^ and/or 
potential NRQCD (pNRQCD)[16] and/or lattice QCD etc at the accuracy for the present 
purposes, whereas, there is no reliable way to offer such accurate values for the color-octet 
non-perturbative matrix elements as those for color-singlet ones even for lattice QCD [13]. 
So far, we may only have rough values of these matrix elements relating to the components 
\(cb)8( 1 So)g) , \(cb)*( 3 S^q) ■ • • by means of the velocity power counting rule of the NRQCD 
effective theory [12]. In order to make a reasonable estimate on the contributions from the 
color-octet components to the hadronic production so as to see if the color-octet components 
are visible experimentally or not, we take a similar way as done in Ref . jl7j| . i.e. the color- 
octet matrix elements are taken to be smaller by certain order in v 2 than those obtained by 
relating them to the color-singlet S-wave wavefunctions at the origin |^(0)| 2 . Whereas, the 
first derivative of the wavefunctions (P-wave) at the origin for the color-singlet components 
with proper factor of m (the reduce mass of the bound state), i.e. 2^151 5 is in a higher 
order in v than the wavefunctions at the origin for the S'-wave functions, "0(0), therefore 
the contributions from the color-octet components \(cb)s( l So)g), \(cb)s( 3 Si)g) ■ ■ ■ (S'-wave) 
to the production of the P-wave excited states hs c and Xb c ma y be comparable with those 
from the color-singlet components |(c6)i( 1 Pi)) and |(cfe)i( 3 Pj)) (P-wave). 

The color-singlet contributions to the P-wave B c meson hadronic production have been 



studied in Refs. 



10 



14. 



18,ll9|. In Ref. 11811 . the hadronic productions of /i R„ a nd y£ were 

nnn 

estimated by applying the fragmentation approach. While in Refs. |10J, HJ, LLi| the authors 
adopted the lowest order (here af) complete approach but only color-singlet components 
|(c6)i( 1 Pi)) and |(c6)i( 3 Pj)) were taken into account, and found that the total cross-section is 
much larger than what predicted by the fragmentation approach [18], and the fragmentation 
results tend to be comparable with those from the a* complete calculation only in the cases 
when the transverse momentum p t of the meson B c is so high as p t > 30 GeV (higher than 



.S-wave production case) 



14, 3]- For the a\ complete estimation, to be complete and to see 



the characteristics, we compute the contributions from the color-octet components to the 
P-wave P c -meson production precisely in the paper. 

The paper is organized as follows: In Section II we will present the basic formulae for 
the lowest order (af) which are used in the complete calculation of the contributions to the 
P-wave hadropro duct ion from the color octet components | (c6) 8 ( 1 5'o)5 f ) an d | (c6) 8 ( 3 S'i)5f) . 
In Section III, we will explain the options about the input parameters, which are necessary 
for the calculations and will present the numerical results properly. In the last section, a 
brief summary and some discussions will be given. 



II. FORMULATION AND TECHNIQUE 

In hadron-hadron collisions at high energy, the gluon-gluon fusion mechanism is the 
dominant one over the others. Hence, to study the hadronic production of the P-wave B c 
mesons, as argued in the preceding section, we need to consider two types of processes, i.e., 
gg — > (c6)i + b + c where (cb)i is in color-singlet 1 P\ and 3 Pj states; and gg — > (c&) 8 + b + c 
where (c6) 8 is in color-octet 1 Sq and 3 Si states. In Ref . . the color-singlet contributions 
to the hadronic production of /i# c and Xb c have been computed, and the useful formulas 
and techniques have been explained there. For simplicity, we will not repeat it here. In the 
present work, our main concern is to complete the calculations on the hadronic production 
of the P-wave B c physical states, especially to pay attention to the contributions from the 
color-octets \(cb)s( l So)g) and |(cfr) 8 ( 3 S'i)g). 

The amplitude for the gg — > (c6) 8 + b + c can be analytically obtained by applying the 
helicity method Q, and (c&) 8 will be 'hadronized' to the state \{cb)s( 1 So)g) or \(cb)s{ 3 Si)g) 
finally. The details of the helicity technique can be found in Ref.|9| and references therein. 
The difference of the present calculations is only in the color configuration from those for 
the leading order B C (B*) meson hadronic production. 

In the case of the hadronic production gg — > (c6)i + b + c where (c6)i is in color-singlet, 
|(c6)i( 1 S'o)) or |(c&)j ( 3 >S'i\) 1 there are only three independent color factors (i.e. three inde- 



pendent color flows 



14 



20]). While for the case of the hadronic production of the color-octet 
states, gg — > (c6) 8 + 6 + c where (c6) 8 is in color-octet and will be hadronized to |(c6) 8 ( 1 S'o)5') 



or | (c6) 8 ( 3 t S , i)5 r ) finally, there are totally ten independent color factors, Cuj (k — 1, 2, • • - , 10), 



where i, j — 1,2, 3 are color indices of the quarks c and b respectively. They are 
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y/2N c ^ >n 

where the indices a and b are color indices for gluon-1 and gluon-2 respectively, and \/2T d 
stands for the color of the color-octet state (cfo)g in the production. N c = 3, for QCD. All 
the color factors in the subprocess can be expressed by the linear combination of the above 
ten independent color factors, i.e., in the amplitude all the color factors may be written in 
terms of these ten explicitly. Thus the total helicity amplitude can be generically expressed 
as 

10 

M^ X ^' X6 \q bl ,q b2 ,q cl ,q c2 ,k h k 2 ) = £ C mij M^ MM > X «\q hU q b2 ,q cU q c2 ik 1 ,k 2 ) , (4) 

m=l 

where Ai, A4, A5 and A6 denote the helicities of the out-going 6-quark, out-going c-quark, 
gluon-1 and gluon-2, respectively, ki and k 2 are the momenta of the gluons; q c i, q b i are the 
momenta of c and b quarks and q c2 , qbi are the momenta of c and b anti-quarks, respec- 
tively. The helicity amplitude M^ 1,A4,A5 ' A6 )(g w , q b2 , g cl , q c2 , k\, k 2 ) can be directly read from 
the Eqs.(4-8) in Ref.Q, with only the color-singlet matrix elements (or the wave function 
at the origin) being replaced by the color-octet matrix elements (0\xlT d "i/j c (a^ H aH)'^lT d Xb\0) 
and (0\x b & l T d %l> c (a' H aH)ijlo'' l T d Xb\fy ■ To get the matrix element squared, one needs first to 
deal with the square of the above ten independent color factors, i.e. (C mi j x C*^) with 
m,n = 1, 2, • • • 10. For reference use, whose values are listed in TABflJ 

According to pQCD factorization theorem, the inclusive B c meson hadropro duct ion can 
be formulated as 

da = J2j dx ij dx 2 F l Hi (x 1 ,/j 1 2 F ) x F^ 2 (x 2 , iJ 1 2 F )da ij ^ {cb)x (x 1 ,x 2 , fi 2 F , ij 2 ,Q 2 ) , (5) 

ij 

where F l Hi (x, fi 2 F ) , F ] H2 (x, fi F ) are parton distribution functions (PDFs) of partons i and j in 
hadrons H\, H 2 , respectively, [i 2 is the 'energy scale squared' where renormalization for the 



TABLE I: The square of the ten independent color factors (including the cross terms), (C m y x C*j ■) 
with m,n= 1, 2, • • • , 10 respectively. 
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subprocess is made; Q 2 is the 'characteristic energy scale of the subprocess squared'; and n 2 F 
is the 'energy scale squared' where the factorization of the PDFs and the hard subprocess 
is made. Usually, in the LO calculation, to avoid the large logarithms the scales of typical 
scale of hard interaction, factorization and 'renormalization' are set to be the same, i.e. 
/i 2 = & = Q 2 . 



III. NUMERICAL RESULTS 



For convenience and as a reference, in the numerical calculation we take the values of the 
radial wave function at the origin and the first derivative of the radial wave function at the 



origin as those given by Refs. 



15 



2l|, i.e., |i?(0)| 2 = 1.54Ge1/ 3 and \R\Q)\ 2 = 0.201GeV 5 
(namely the two values roughly mean v 2 ~ 0.1, if the reduced mass m ~ 1.2GeV), which 
relate to the non-perturbative matrix elements of the color singlet production definitely. 
Whereas there is no reliable way to compute the production color-octet matrix elements 
(0|xiT d ^ c (a^a^)^T d Xfe|0) and (Olx^T^^a^a^)^^^^^). Although we do not know 
the exact values of the above two octet matrix elements, according to NRQCD scale rule we 



know that they are at the same order of the P-wave color-singlet matrix elements, which are 
one order in v 2 higher than the iS-wave color-singlet matrix elements accordingly. Moreover 
according to the heavy-quark spin symmetry 12], the values of the above two production 
matrix elements have the approximate relation: 

(0\xWT d MaUH)^T d Xb \0) = (2J+ 1) • {^ h T d ^ H a H )^ c T d Xb\m + 0(v% (6) 
where J is the total angular momentum of the hadron state. In our numerical calculations, 



as done in Ref. 17|, the values of the color-octet matrix elements are set to be smaller by 
certain order v 2 than those obtained by relating to the S-wave wave functions at the origin 
^(O)! 2 for the color singlet. More specifically, based on the velocity scale rule, we can 
estimate 



(Q\xl^\B c CS )) . [l + 0(v 4 )} , (7) 



where the second equation comes from the vacuum-saturation approximation. As(v) is of 
the order v 2 or so, and we take it to be within the region of 0.10 to 0.30, which is in consistent 
with the identification: As(v) ~ a s (Mv) and has covered the possible variation due to the 
different ways to obtain the wave functions at the origin (S-wave) and the first derivative of 
the wave functions at the origin (P-wave) etc. 

As argued in Ref. the equations: P = q c \ + g&2, Qci = m c> 1b2 = m b an d P 2 = M 2 
must be satisfied simultaneously, and furthermore, only when these equations are satisfied 
can the gauge invariance of the amplitude be guaranteed. Hence, in this paper we proceed the 
numerical calculation by making the choice: the masses of c and b quarks are m c = 1.50 GeV 
and m b = 4.90 GeV, respectively; the mass of the bound states |(c&)i( 1 P 1 )), |(c6)!( 3 Pj)), 
|(c6)s( 1 5'o)^) and | (cfyg^S^g) are therefore to be M = 6.40 GeV, the sum of the two heavy 
quark masses. 

There are a couple of uncertainty sources in the theoretical estimates on the B c meson 
hadronic production, such as those from the strong coupling constant, the oe s ; from the PDFs 
of different sets of fittings etc.. In the present calculations, the factorization energy scale is 
set to be the transverse mass squared of the (cb) bound states, i.e., Q 2 = M 2 = [M 2 + p 2 ) 
where pt is the transverse momentum of the bound state; the PDF of version CTEQ6L |2jJ 
and the leading order a s running above Aq^ 4 ^ = 0.326 GeV are adopted. 



TABLE II: Total cross-section (in unit of nb) for the hadronic production of the (cb) meson at 
LHC (14.0 TeV) and TEVATRON (1.96 TeV), where for short the \?S )i) denotes (cb) state in 
color-singlet ( 1 5'o) configuration, and so forth. Here mj = 4.90 GeV, m c = 1.50 GeV and M = 6.40 
GeV. For the color-octet matrix elements, we take As(v) G (0.10, 0.30). 





\Cs )i) 


\esih) 


\( 1 S )sg) 


\( 3 Si)s 9 ) 


\ePi)i) 


K 3 ^)!) 


l( 3 ^i)i) 


ICftJi) 


LHC 


71.1 


177. 


(0.357, 3.21) 


(1.58, 14.2) 


9.12 


3.29 


7.38 


20.4 


TEVATRON 


5.50 


13.4 


(0.0284, 0.256) 


(0.129, 1.16) 


0.655 


0.256 


0.560 


1.35 




P t (GeV) y 

FIG. 1: Distributions in pt and y of the hadronic production of (cb) meson at LHC. The dashed line, 
solid line, dash-dot line, dotted line represent the color-singlet Pi, 3 Po, 3 -Pi, and 3 P2, respectively. 
The lower and upper shaded bands stand for the color-octet 1 Sq and 3 Si states respectively, whose 
upper limit corresponds to As(v) = 0.3 and lower limit corresponds to A$(v) = 0.1. 

In TABLE|nJ we show the total cross-sections for hadronic production of the (cb) meson 
at LHC and TEVATRON respectively, where the (cb) may mean in |(c6)i( 1 Pi)), \(cb)i( 3 Pj)) , 
I (cb)s( 1 So)g) and \(cb)g( 3 Si)g) configurations respectively, and the results of the |(c6)i( 1 S'o)) 
(the dominant component for B c ) and |(cfo)i( 3 5'i)) (the dominant component for B*) pro- 
duction are included for comparison. In the table, we take A s (v) G (0.10, 0.30). One may 
observe that the cross-sections of the color-octet S-wave states are comparable to those of 
the color-singlet P-wave states, and they are in the same order in v 2 expansion as expected 
by NRQCD. 




FIG. 2: Distributions in pt and y of the hadronic production of (cb) meson at TEVATRON. 
The dashed line, solid line, dash-dot line, and dotted line represent the color-singlet Pi, 3 Po, 
3 Pi and 3 P2; respectively. The lower and upper shaded bands stand for color-octet 1 Sq and 3 S± 
states respectively, whose upper limit corresponds to As(v) = 0.3 and lower limit corresponds to 
A s (v) = 0.1. 

To see this point more clearly, we calculate the distributions of the corresponding bound 
states on the transverse momentum p t and the rapidity y respectively, and plot the results 
in Figs. (jlJ2)) . From the figures Figs.(^Q EJ, one may observe that the contributions to the 
total cross section from the color-octet component are in comparable with those from the 
color-singlet component in the P-wave B c meson hadronic production. However, in the 
Pt distribution the color-octet contributions drop more rapidly than those from the color- 
singlet parts, which exhibits more clearly through looking at the ratio, Riptcut) — a " ctet , 

^singlet 

where a oc t e t stands for the total contribution from the color-octet and <J S i n giet for that from the 
color-singlet. The variation of the dependence of R(ptcut) on p tcu t is presented in TABLE llTTl 
Here, for simplicity, we take a middle value of A|(f), i.e. A 2 s {v) = 0.05 (or A s (v) = 0.224), 
for the color-octet matrix elements. 

Since there is no detector in high energy hadronic collisions which can directly detect all 
the events, especially those with a small p t and/or a large rapidity y, so for experimental 
studies and for practical purpose, only the events with proper kinematic cuts on p t and/or y 
are taken into account. Now let us compute the observables with cuts on p t and/or y, so as 
to see the the cut effects in h Bc and Xb c production precisely. Namely we have investigated 
the effects in production rate and various differential cross sections and for the color-octet 



TABLE III: The correlation of R(p 

tout) — °" octet with the Ptcut- F° r color-octet matrix elements, 

^singlet 

we take Ag(t)) = 0.05, which is in the middle of the range (0.01,0.09). 
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FIG. 3: pt distributions of various contributions in hs c and Xb c production, with several different 
cuts in rapidity (y C ut) at LHC (left) and at TEVATRON (right) energies. Dashed line (next to top) 
with y cu t = 2.0; dash-dot line (middle) with y cu t = 1.5; dotted line (next bottom) with y cu t = 1.0; 
diamond line (bottom) with y cut = 0.5 and solid line (top) without y cut . The color-octet matrix 
elements are set at A 2 s (v) = 0.05. 

production and the color-singlet production of S"-wave (B c and B*) and P-wave (hs c and 
Xb c ) respectively. Explicitly, we take A|(i>) = 0.05 for the color-octet matrix elements in 
the investigation. Considering the abilities on measuring p t and y rapidity of B c for CDF, 
DO, BTeV at TEVATRON and for ATLAS, CMS, and LHC-B at the LHC, we compute 
the p t distributions with the rapidity cuts y cut = 1.5, and the y distributions with the 
transverse momentum cut p tcu t = 5 GeV accordingly. The results with four rapidity cuts, 
y cu t = (0.5,1.0,1.5,2.0), are put together in Fig.® and the results with four transverse 
momentum cuts, p tcut = 5, 20, 35, 50 GeV, are put in Fig.®. 



FIG. 4: y distributions of all the considered contributions, i.e., from both color-octet S'-wave states 
and color-singlet P-wave states, with various transverse momentum cuts (y C ut) at LHC (left) and 
at TEVATRON (right) energies. Here, the Diamond line represents what with ptcut = 5 GeV; 
dash-dot line with ptcut = 20 GeV; dotted line with ptcut = 35 GeV; dashed line with ptcut = 50 
GeV, and the solid line with ptcut = 0. Of the color-octet matrix elements, we take A|(f ) = 0.05. 

IV. DISCUSSIONS AND SUMMARY 

In the paper, according to the NRQCD expectation about the importance for the compo- 
nents, we have precisely investigated the contributions of the S-wave color-octet components 
to the hadronic production of the P-wave B c states. Our final results show that contribu- 
tions from the color-octet ones are comparable to those from the color-singlet ones in the 
production of the P-wave B c excited states if the scale rule of NRQCD works well and the 
value of the relative velocity squared v 2 is really in the possible region 0.1 ~ 0.3, as indicated 
by potential models for the (cb) binding system for instance. Therefore, to make a soundly 
prediction to the hadronic production of the P-wave (cb) meson at leading order in v 2 , one 
needs to take both contributions from the color-octet S'-wave Fock states and those from 
the color-singlet P-wave Fock states into account. With A|(w) = 0.05 and ptcut = 0.0 GeV, 
the total contributions from the color-octet S'-wave Fock states to the total production cross 
section is about ~ 22% from the color-singlet P-wave states at LHC and about ~ 25% at 
TEVATRON respectively, and such contributions decrease with the increment in ptcut, e.g. 
at Ptcut = 50 GeV, it reduces to ~ 4% at both energies of LHC and TEVATRON. This 
character may used for people to distinguish the color-singlet contributions from those from 



color-octet ones in future when there are enough hs c and Xb c data. 

Note that for shortening the paper and not changing the main feature, here we have not 
taken the possible mixing between hs c and Xb =1 account at all. According to NRQCD 
scale rule, except the components \ (cb)g( 1 So)g)) and Kcfys^S^g)), there is no enhancement 
factor, such as that comes from the wave function at the origin for the S-wave components 
virus the derivative of the wave function at the origin for the leading P-wave color-singlet 
ones, so we consider the contributions from the other high order Fock components in Eq. (j2J) 
small, and ignore them at all. 

In summary, the total cross-section of the h Bc and Xb c production, including both singlet 
and octet contributions, at the lowest order in a s and v 2 expansions, may be so bigger than 
a half of the direct production rate of the ground state B^Sq) (see TABLE. II, roughly 
speaking: ~ 0.8 for LHC and ~ 0.7 for TEVATRON). Considering the fact that almost 
all of the low-laying excited states decay to the ground state B c , it also means that the 
full low-laying P-wave state production 'promptly' contributes the P c ( 1 S'o) production by a 
bigger factor than 0.5 of the direct one. Especially, if A|(t>) ~ 0.05, the contributions from 
color-octet components themselves may contribute promptly a factor ~ 0.1 of the direct 
hadronic production of B c meson. Furthermore, our complete calculations also indicate 
that the observation of the P-wave B c states is not very far from the present experimental 
capability in principle. 
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